Differentiating dark energy and modified gravity with galaxy redshift surveys 
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The observed cosmic acceleration today could be due to an unknown energy component (dark 
energy), or a modification to general relativity (modified gravity). If dark energy models and 
modified gravity models are required to predict the same cosmic expansion history H(z), they will 
predict different growth rate for cosmic large scale structure, f g (z). If gravity is not modified, the 
measured H(z) leads to a unique prediction for f g (z), f g (z), if dark energy and dark matter are 
separate. Comparing f^(z) with the measured f g (z) provides a transparent and straightforward 
test of gravity. We show that a simple x 2 test provides a general figure-of-merit for our ability 
to distinguish between dark energy and modified gravity given the measured H(z) and f g (z). We 
find that a magnitude-limited NIR galaxy redshift survey covering >10,000 (deg) 2 and the redshift 
range of 0.5 < z < 2 can be used to measure H(z) to 1-2% accuracy via baryon acoustic oscillation 
measurements, and f g (z) to the accuracy of a few percent via the measurement of redshift-space 
distortions and the bias factor which describes how light traces mass. We show that if the H(z) 
data are fit by both a DGP gravity model and an equivalent dark energy model that predict the 
same H{z), a survey area of 11,931 (deg) 2 is required to rule out the DGP gravity model at the 
99.99% confidence level. It is feasible for such a galaxy redshift survey to be carried out by the next 
generation space missions from NASA and ESA, and it will revolutionize our understanding of the 
universe by differentiating between dark energy and modified gravity. 

3.5 



I. INTRODUCTION 

The observed cosmic acceleration today 0, H| could 
be due to an unknown energy component (dark energy, 
e -9-i 0)i or a modification to general relativity (modi- 
fied gravity, e.g., [J, Q). Ref.pj] contains reviews with 
more complete lists of references. Illuminating the na- 
ture of dark energy is one of the most exciting challenges 
in cosmology today. 

The cosmic expansion history, H(z) = (d\na/dt) (a 
is the cosmic scale factor), and the growth rate for 
cosmic large scale structure, f g (z) = dhiS/dliia [5 — 
(Pm — Pm) I Pro] ■ are t w0 functions of redshift z that can 
be measured from cosmological data. They provide in- 
dependent and complementary probes of the nature of 
the observed cosmic acceleration 0, Q. The precisely 
measured H{z) and ft m lead to a unique prediction for 
f g {z) in the absence of modified gravity, f g {z), if dark 
energy and dark matter are separate. Comparing f^{z) 
with the measured f g {z) obs provides a transparent and 
straightforward test of gravity (see Fig.l). If gravity is 
not modified, H(z) and f g (z) together provide stronger 
constraints on dark energy models [loj . 

Using the VVDS data, Ref. [ill] demonstrated that a 
magnitude-limited galaxy redshift survey can be used to 
measure f g (z) via measurements of redshift-space distor- 
tion parameter 
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FIG. 1: Current and expected future measurements of the cosmic 
expansion history H(z) = HqE(z) and the growth rate of cosmic 
large scale structure f g (z) = dlnS/dlna (S = (p m — pm)/pm)i 0, 
is the cosmic scale factor). Note that the fiducial model assumed 
for the future galaxy redshift survey is a dark energy model with 
the same H (z) as that of the DGP model. These two models have 
identical expansion histories H{z) [solid line in panel (a)], but very 
different growth rates f g (z) [solid and dashed lines in panel (b)]. 



and the bias parameter b(z) (which describes how light 
traces mass) from galaxy clustering. In this paper 
we show that a feasible, sufficiently wide and deep 
magnitude-limited galaxy redshift survey will allow us 



to unambiguously differentiate between dark energy and 
modified gravity by providing precise measurements of 
H(z) and f g (z) (see Fig.l). 
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II. MODELS 



If the present cosmic acceleration is caused by dark 
energy, E(z) = H(z)/H = [fl m (l + z) 3 + fi fe (l + zf + 
QxX(z)] 1 / 2 , where X(z) = p x (z)/p x (0), with p x (z) 
denoting the dark energy density. The linear growth rate 
f g = dhiDi/dhia is determined by solving the equation 
for D x = 5^(x,t)/5(x), 



D>;(t) + 2E(z)D[(t) 



^fl m (l + z) 3 D 1 = 0, (2) 



H z - - = 



H 

'Y) 
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(3) 



where ro is a parameter in DGP gravity, and p, 
p m (0)(l + zf. Solving Eq.® gives 



Af; 



where primes denote d/d(Hot), and we have assumed 
that dark energy and dark matter are separate. 

In the simplest alternatives to dark energy, the present 
cosmic acceleration is caused by a modification to general 
relativity. The only rigorously worked example is the 
DGP gravity model [EI3]j which can be described by a 
modified Friedmann equation: 
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with fi^ ee p m (0)/p% p° c = 3H$/(8irG). The added 
superscript "0" in denotes that this is the matter 
fraction today in the DGP gravity model. Note that 
consistency at z = 0, -E'(O) = 1 requires that 
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FIG. 2: The accuracy of approximate expressions for f g (z) for 
various models. 



imply that Q 



eff 
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so the DGP gravity model is parametrized by a single 
parameter, f2^. The linear growth factor in the DGP 
gravity model is given by Q 



D'{{T)+2E(z)D[{T)~-n m (l+z) 3 D 1 ( 1 
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where 



O-DGP 



1 - 2H r + 2(g r ) 2 
1 - 2H r 



= 0, 
(6) 

(7) 



The dark energy model equivalent of the DGP gravity 
model is specified by requiring 



8nGp% f H 
3 r ' 



(8) 



Eq.® and the conservation of energy and momentum 
equation, 



p e /J+Hp e /J+P e /e f )H = Q, 



(9) 



where 



_ 8nGp m (z) _ Q° m {l + zf 
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(11) 

As a — ► 1, 
This means 



that the matter transfer function for the dark energy 
model equivalent of viable DGP gravity model (fi^, < 0.3 
and w < —0.5) is very close to that of the ACDM model 
at k > 0.001 hMpc' 1 . [12 

It is very easy and straightforward to integrate Eqs.([2|) 
and ^ to obtain f g for dark energy models and DGP 
gravity models, with the initial condition that for a — > 
0, Di(a) = a (which assumes that the dark energy 
or modified gravity are negligible at sufficiently early 
times). There are well known approximations to f g , 
with f g (z) = Q m (a)^ n for dark energy models [H[ , and 
f g (z) = n m (af/ 3 for DGP gravity models 0. Fig.2 
shows that these powerlaw approximations of f g are 
not sufficiently accurate for future galaxy redshift sur- 
veys that can measure f g to a few percent accuracy in 
Az = 0.2 redshift bins. 
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III. ANALYSIS TECHNIQUE 

Galaxy redshift surveys allow us to measure both H{z) 
and f g (z) through baryon acoustic oscillation (BAO) 
measurements [l5 , [l6 jl7l . l32j and redshift-space distor- 
tion measurements [1 l| . BAO in the observed galaxy 
power spectrum have the characteristic scale determined 
by the comoving sound horizon at recombination, which 
is precisely measured by the cosmic microwave back- 
ground (CMB) anisotropy data Comparing the ob- 
served BAO scales with the expected values gives H(z) 
in the radial direction, and Da(z) [the angular diameter 
distance Da(z) = r(z)/(l + z), where r(z) is the coor- 
dinate or comoving distance] in the transverse direction. 
We will only estimate the accuracy to which H(z) and 
fg(z) can be determined from galaxy redshift surveys in 
dark energy models (the error bars in Fig.l). 1 

The observed power spectrum is reconstructed using 
a particular reference cosmology, including the effects of 
bias and redshift-space distortions [16| : 



Pobs > ^) 



[ DA (zrff h(z) b 2 _ 

[D A (z)] 2 H(zY-f V 
G(z)~\ 2 

~qJq^ Pmatter{k)z=0 + Pshoti^) 



where the growth factor G(z) and the growth rate 
f g (z) = (3b(z) are related via f g (z) = d\nG(z)/d\na, 
and fi = k • r/k, with f denoting the unit vector along 
the line of sight; k is the wavevector with |k| = k. Hence 
fi 2 = k 2 /k 2 = + k 2 ). The values in the refer- 

ence cosmology are denoted by the subscript "ref " , while 
those in the true cosmology have no subscript. Note that 



<„ re f 



Da{z ) 
D a (z) 



ref ■ 
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II 



H(z) 



ref 



H(z) 



(13) 



Eq. (|12]) characterizes the dependence of the observed 
galaxy power spectrum on H(z) and Da(z) due to BAO, 
as well as the sensitivity of a galaxy redshift survey to 
the redshift-space distortion parameter (3 |2l| . 

To study the expected impact of future galaxy redshift 
surveys, we use the Fisher matrix formalism. In the limit 
where the length scale corresponding to the survey vol- 
ume is much larger than the scale of any features in P(k), 
we can assume that the likelihood function for the band 
powers of a galaxy redshift survey is Gaussian [221 ] . Then 
the Fisher matrix can be approximated as [23[ 



<91nP(k) SlnP(k) 



dpi 



d Pj 



Veff(k) 



2(2tt) 3 



(14) 



where pi are the parameters to be estimated from data, 
and the derivatives are evaluated at parameter values of 
the fiducial model. The effective volume of the survey 



V e ff(k,n) 



dr A 



n(r)P(fc, fi) 
n(r)P(fc,/i) + 1 



nP(fc, /i) 
nP(k,fi) + 1 



V, 



survey 3 



(15) 



where the comoving number density n is assumed to 
only depend on the redshift for simplicity. Note that 
the Fisher matrix Py is the inverse of the covariancc 
matrix of the parameters pi if the pt are Gaussian dis- 
tributed. Eq. (|14p propagates the measurement error in 
lnP(k) (which is proportional to [^//(k)] -1 / 2 ) into mea- 
surement errors for the parameters Pi . 

Since we do not include nonlinear effects, we only con- 
sider wavenumbers smaller than a minimum value of 
non-linearity. Following [l5j . we take k m i n — 0, and 
k max given by requiring that the variance of matter fluc- 
tuations in a sphere of radius R, tr 2 (P) = 0.35, for 
R = 7r/ (2k max ). We will also give results for <r 2 {R) = 0.2 
for comparison. In addition, we impose a uniform upper 
limit of k max < 0.2/iMpc" 1 , to ensure that we are only 
considering the conservative linear regime essentially un- 
affected by nonlinear effects. [25| shows that nonlinear 
effects can be accurately taken into account. 26] shows 
that the BAO signal is boosted when these effects are 
properly included in the Hubble Volume simulation. We 
assume fl b = 0.045, h = 0.7, 6=1, and nP = 3 [l|; 
this is conservative since nP > 3 at any redshift for a 
magnitude-limited survey. 

The observed galaxy power spectrum in a given red- 
shift shell centered at redshift zi can be described by a set 
of parameters, {H(zi), D A (z l ) 1 G(zi), /3(z,), P l shoV n s , 
u> m , Wf,}, where rig is the power-law index of the primor- 
dial matter power spectrum, u> m = fl m h 2 , and uji, = fl h 2 
(h is the dimensionless Hubble constant). Note that P(k) 
does not depend on h if k is in units of Mpc -1 , since the 
matter transfer function T{k) only depends on u) m and 
[13] > 2 ^ the dark energy dependence of T(k) can be 
neglected. Since G(z), b, and the power spectrum nor- 
malization Pq are completely degenerate in Eq. (|12|) . we 



have defined G{ Zi ) = bG{z) P 1/2 /G(0). 

The square roots of diagonal elements of the inverse 
of the full Fisher matrix of Eq. (fT^|) gives the estimated 
smallest possible measurement errors on the assumed pa- 
rameters. The parameters of interest are {H(zi), Da{zi), 
/3(zi)}, all other parameters are marginalized over. Note 
that the estimated errors we obtain here are independent 



1 gives a more precise treatment of redshift-space distortions, 
and |19ll studies power spectra in alternative gravity models. 



2 Massive neutrinos can suppress the galaxy power spectrum am- 
plitudes by > 4% on BAO scales |24| . It will be important for 
future work to quantify the effect of massive neutrinos on the 
measurement of H(z) and f g (z). 
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of cosmological priors, 3 thus scale with (area)^ 1 / 2 for a 
fixed survey depth. 

Fig.l shows the errors on H(z) and f g (z) = B{z)b{z) 
for a dark energy model that gives the same H(z) as a 
DGP gravity model with the same f2^n, for a rcdshift 
survey covering 11,931 (deg) 2 , and the redshift range 
0.5 < z < 2 [<7 2 {R) = 0.35 assumed]. Note that the 
Da{z) measured from the same redshift survey provides 
additional constraints on H{z) that can be used for cross- 
checking to eliminate systematic effects. We have ne- 
glected the very weak dependence of the transfer function 
on dark energy at very large scales in this model [l2l ]. 
and added Aln6 = 0.01 {(area)/[28,600 (deg) 2 ]}" 1 / 2 in 
quadrature to the estimated error on (3. A 

[36| developed the method for measurin g b( z) from the 
galaxy bispectrum, which was applied by [37] to the 2dF 
data. Assuming that [HI 



S g = M(x) + i6 2 5 2 (x), 



(16) 



the galaxy bispectrum 



<^kAk 2 <W = (2 7 r) 3 |p 9 (k 1 )P 9 (k 2 ) 



J(ki,k 2 ) 



h 
b 2 



eye. U^Ckx + ka + ka), 



(17) 



where J is a function that depends on the shape of the 
triangle formed by (ki, k 2 , in k space, but only de- 
pends very weakly on cosmology 36] . 

Ref. [HI used Monte Carlo N-body simulation to study 
the extraction of the BAO scales. For comparison, we 
calculated {H(zi), D^fe)} for the same fiducial model 
as considered by [l5j (with the same assumptions and 
cutoffs in k) , and obtained results that are within 30% of 
the values given by the fitting formulae from |28j| . This 
is reassuring, as it validates the approach of using the 
Fisher matrix formalism to forecast the parameter accu- 
racies for future redshift surveys. 5 



IV. OBSERVATIONAL METHODS 

H(z) can be probed using multiple techniques. It can 
be measured using Type la supernovae (SNe la) as cos- 
mological standard candles [29| . CMB and large scale 



3 Priors on <j m , ui^, fife, and ng will be required to obtain the 
errors on dark energy parameters. 

4 This A In b estimate comes from extrapolating 2dF measurement 
of b = 1.04 ± 0.11 at z ~ 0.15 for an effective survey area of 
1300x127000/245591=672 (deg) 2 [37|], and assuming a factor of 
1.6 improvement for a NIR space mission that can detect galax- 
ies at a much higher number density. This A In b estimate is 
comparable (and larger) than that estimated by |38| for imaging 
surveys at z < 2. 

5 Ref. [23 found similar agreement in their comparison. 



structure data provide constraints on cosmological pa- 
rameters that help tighten the constraints on H(z) [30]. 
Fig.l(a) shows the H(z) given by Eq.Q with = 0.25 
(solid line) , as well as a cosmological constant model with 
O m = 0.3, Qa = 0.7 (dotted line). Clearly, both these fit 
the constraints on H (z) from current data [30] (no priors 
assumed). 6 

BAO measurements from a very wide and deep galaxy 
redshift survey provide a direct precise measurement of 
H(z) [see Fig. 1(a)]. Suppose H(z) is measured to be 
H 2 — H/r = 87rGp m /3 [see solid line in Fig. 1(a)] and fi m 
is known accurately, Eq. ^ yields a unique prediction for 
f g {z), f^{z), assuming that gravity is not modified [see 
the dashed line in Fig. 1(b)]. 

The measurement of f g [z) can be obtained through 
independent measurements of 8 = f g (z)/b and b(z) [111 ] . 
The parameter [3 can be measured directly from galaxy 
redshift survey data by studying the observed redshift- 
space correlation function [33l. l34j . The bias factor b(z) 
can be measured by studying galaxy clustering properties 
(for example, the galaxy bispectrum) from the galaxy 
redshift survey data [37[. Independent measurements of 
Biz) and b{z) have only been published for the 2dF data 

Fig. 1(b) shows the f g (z) for the DGP gravity model 
with Q® n — 0.25 (solid line), as well as a dark energy 
model that gives the same H(z) for the same (dashed 
line). The cosmological constant model from Fig. 1(a) is 
also shown (dotted line). Clearly, current data can not 
differentiate between dark energy and modified gravity. 

A very wide and deep galaxy redshift survey provides 
measurement of f g (z) accurate to a few percent [see 
Fig. 1(b)]; this will allow an unambiguous distinction be- 
tween dark energy models and modified gravity models 
that give identical H (z) [see the solid and dashed lines in 
Fig. 1(b)]. A simple x 2 test can provide a general figure- 
of-merit for our ability to distinguish between dark en- 
ergy and modified gravity models that fit the measured 
H(z) but predict different f g (z). If the measurement 
errors are normally distributed, A\ 2 = X 2 ( s ) — x 2 ( s o) 
is distributed as a chi-square distribution with n de- 
grees of freedom (n is the number of data points) , where 
s is the test model, and So is the bestfit model mea- 
sured from data. P(x 2 \ n ) = 99.99% corresponds to 
Ax 2 = 29.877 for n = 7. Assuming that x 2 ( s o) = n, 
we find that x 2 ( s ) — 36.877. In Fig.l, we assume that 
the true model is a dark energy model with fi^ = 0.25, 
H 2 - H/r = 87rGp m /3, with Hr = 1/(1 - n°J. For 
a linear cutoff given by <J 2 (R) = 0.35 (or 0.2), a survey 
covering 11,931 (deg) 2 would rule out the DGP gravity 
model that gives the same H{z) and at 99.99% (or 
95%) C.L.; a survey covering 13,912 (deg) 2 would rule 
out the DGP gravity model at 99.999% (or 99%) C.L.. 



Ref. [30j uses WMAP three year data [3], 182 type la supernovae 
[3l| . and the SDSS baryon acoustic scale measurement 
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V. CONCLUSIONS 

Discovering the nature of dark energy has been iden- 
tified as a high priority by both NASA and ESA. A 
magnitude-limited NIR galaxy redshift survey, covering 
>10,000 (deg) 2 and the redshift range 0.5 < z < 2, can be 
feasibly carried out by a space mission that uses MEMS 
technology to obtain 5000-10,000 galaxy spectra simulta- 
neously [4(| IH| • The low background from space enables 
very short exposure times to obtain galaxy spectra to 
z ~ 2, making it practical to carry out a magnitude- 
limited NIR galaxy redshift survey over > 10,000 (deg) 2 
in only a few years. A magnitude-limited galaxy red- 
shift survey over > 10,000 (deg) 2 will enable robust and 
precise determination of b(z) using multiple techniques 
and with sufficient statistics [ill HE H3 ■ This is criti- 
cal for determining f g (z) using measurements of redshift- 
space distortions. Such a survey will also enable rigorous 
study of the systematic uncertainties of BAO, and accu- 
rate measurements of redshift-space distortions. 

Ref.0 studied the use of weak lensing shear maps 
to differentiate between dark energy and modified grav- 
ity, complementary to what we have studied in this pa- 
per. While both weak lensing surveys and galaxy red- 
shift surveys can provide accurate measurements of H(z) 
(if the systematic uncertainties are properly modeled 

' Such a survey would allow us to distinguish between dark energy 
and modified gravity even if dark energy is clustered such that 
f g , bias, and redshift distortions are scale-dependent |42| , since a 
dark energy model and a modified gravity model generally have 



and controlled), galaxy redshift surveys can potentially 
provide the most accurate measurement of f g (z) [com- 
pare Fig. 2 of H with Fig.l of this paper, noting that 
f g (z) = dlnG(z)/dlna\. 

We have shown that a magnitude-limited NIR galaxy 
redshift survey covering >10,000 (deg) 2 and 0.5 < z < 2 
can provide precise measurements of the cosmic expan- 
sion history H(z), and the growth rate of cosmic large 
scale structure f g (z). These provide model-independent 
constraints on dark energy and the nature of gravity. The 
precisely measured H(z) can be used to predict f g {z) 
expected in the absence of modified gravity, f^[z), if 
dark energy and dark matter are separate. Comparing 
fg{z) with f g (z) obs provides a transparent and pow- 
erful probe of modified gravity. This will allow us to 
illuminate the nature of the observed cosmic accelera- 
tion by differentiating between dark energy and modified 
gravity [see Fig.l]. A magnitude-limited survey covering 
11,931 (deg) 2 can rule out the DGP gravity model at the 
99.99% confidence level. 7 If this technologically feasible 
survey is carried out by a space mission, it will have a 
revolutionary effect on our understanding of the universe. 

Acknowledgements: I thank Gigi Guzzo for helpful 
comments on a draft of this paper, and Chris Blake, Craig 
Wheeler, and Eiichiro Komatsu for useful discussions. 

different redshift dependences of the modified growth rate, and 
the data of such a survey can be analyzed in multi redshift slices, 
on multi scales, and using different populations of galaxies. 
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